Dirac-Weyl electrons in a periodic spin-orbit potential 
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Graphene super-structures have been widely studied but the original form of the SU (2) Hamilto- 
nian was never modified. We study SU(2)(8)SU(2) super-structures arising from spatial modulation 
of spin-orbit interactions and derive an analytic band condition valid for a lattice momentum along 
the direction of modulation of the spin-orbit interactions. The simple form of this band condition 
enables us to estimate the size of gaps due to avoided band crossings and gives insight into the 
dependence of the band structure on the width of the potential. We also investigate band structures 
for the case where the lattice momentum forms a finite angle (j) with respect to the modulation 
direction of the spin-orbit interactions. 

PACS numbers: 61.48.De,72.80.Vp,72.25.Mk 



A(x) X(x) 

^Ld □ □ n n n , 



X 




a b i 



FIG. 1. Sketch of the graphene plane in the presence of a 
modulation of the SOIs along the x-direction. The periodicity 
length is ^ = a + 5 with a the length of the region with SOI. 



I. INTRODUCTION 

Graphene, a two- dimensional (2D) network of s'p^ car- 
bon atoms forming a honeycomb lattice, is well known 
for its low-energy electronic excitations. These can be 
described as Dirac-Weyl quasi-particles, characterized by 
a linear spectrum and a chiral nature leading to many of 
the unusual properties of this new material [1]. In the 
present letter we present yet another unusual effect, aris- 
ing from the special structure of the Rashba Hamiltonian 
for graphene. We impose a super-lattice of periodically 
modulated spin-orbit interactions (SOIs) on graphene. 
As SOIs interact with the spin degree of freedom, the 
original SU(2) structure of the graphene Hamiltonian is 
modified to a SU(2)(g)SU(2) structure. For such super- 
lattices it is usually not possible to obtain a simple band 
equation, unlike in the SU(2) case, where analytic band 
equations were already obtained for a large number of 
different cases [2-6]. Despite this, for a lattice momen- 
tum along the modulation direction of the SOI, we derive 
analytic band equations, which look very similar to those 
obtained for the easier SU (2) case. 

Graphene has two possible mechanisms of SOIs: the 
extrinsic and the intrinsic one. The intrinsic SOI origi- 
nates from carbon intra-atomic SOI. It opens a gap in the 
energy spectrum and converts graphene into a topologi- 
cal insulator with a quantized spin-Hall effect [7]. This 
term has been estimated to be rather weak in clean flat 
graphene [8-11]. However, a recent theoretical proposal 



shows that the presence of indium and thallium adatoms 
can enhance the gap associated to this effect. Even for 
low coverages of 6%, the gap for indium is of the or- 
der of 100 K, while for thallium it approaches room 
temperature [12]. The extrinsic Rashba-like SOI origi- 
nates from interactions with the substrate, presence of a 
perpendicular external electric field, or curvature of the 
graphene membrane [8-10, 13, 14]. A series of recent 
experiments [15-17] have reported the possibility to con- 
trol and vary the Rashba term by opportune substrate 
engineering. A combination of the latter method and a 
nanomasking technique for controlling the adatom depo- 
sition would make variations of both the intrinsic and 
extrinsic SOI strengths feasible on sub-micrometer scales 
and thus our results applicable for experiments. 

Modulated SOI leads to spin-polarized energy bands: 
this is essential for realizing graphene-based spintron- 
ics devices [18] and especially for generating pure spin- 
currents, using, e.g.^ quantum pumping [19] or the spin 
ratchet effect [20-22]. 

In the present work we introduce the spin-dependent 
transfer matrix (TM) method and use it to derive the 
spin-polarized energy bands in the case of SOI modu- 
lation. These simple, analytic band equations allow to 
directly obtain properties characterizing the band struc- 
ture, say the band gaps and the dependence on the length 
of the SOI potential. 

II. MODEL AND FORMALISM 

We use the single valley Dirac-Weyl Hamiltonian [1], 
which in absence of SOIs, reads 

no=VFcr-p, (1) 

where ~ lO^m s~^ is the Fermi velocity, a is the 
vector of Pauli matrices associated with the pseudo-spin 
and p = —\h{dx-,dy) the momentum operator in the 2D 
graphene plane. In the following we use units in which 
fi — vy = 1. 

As mentioned before two types of SOIs are allowed 
in graphene: the intrinsic and the extrinsic or Rashba 
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SOI [7]. The intrinsic SOI Hamiltonian is derived from a 
second nearest neighbor hopping term [8, 13, 14]. In the 
low-energy limit the Hamiltonian reads 



ni = A(x)cr^5^ 



(2) 



where A{x) is a space dependent coupling constant, 
which can be implemented by appropriate deposition of 
adatoms. Here Sz is the Pauli matrix associated with the 
spin of the Dirac-Weyl fermions. The intrinsic SOI does 
not lift spin degeneracy but opens a gap in the energy 
spectrum [7]. 

When inversion symmetry is broken, a Rashba SOI 
term is possible. It reads 



\{x) 



(3) 



Here, \{x) is a space dependent coupling. Contrary to 
1-Li this term lifts spin degeneracy but does not open a 
gap [7, 23]. 

The single valley model (1) is valid within the approx- 
imation that the length scales of variation of the SOIs — 
via the \{x) and A(x) — are assumed to be much larger 
than graphene's lattice constant ao ~ 1.46 A but much 
smaller than the typical Fermi wavelength of quasiparti- 
cles Ap. Since close to the Dirac points Ap ~ l/kl, at low 
energy e this approximation is justified. This assumption 
ensures that the variation of SOIs, close to the K points, 
can be approximated as a sharp change. Therefore we 
choose A(x) and A(x) to be periodic, piece-wise constant 
functions which are either zero or have magnitude A and 
A, respectively (c. f. Fig. 1). 

In the following we introduce a generalization of the 
well known TM method [2]. It is valid when the lattice 
momentum is fixed along the modulation direction and 
leads to analytic band equations. For general directions 
of the lattice momentum, as discussed later on, we will 
consider a different method based on the Bloch theorem. 



III. SPIN-DEPENDENT TM METHOD 

We consider a unit cell of length ^ = a + 6 with a region 
h without SOI and a region a with finite SOI as depicted 
in Fig. 1. In the nth unit cell the wave function 
can be written as a superposition of the forward (F) and 
backward (B) propagating spin up (t) and spin down {\) 
spinor solutions of the stationary Schrodinger equation 

Vt = {0F,t5 0B,t5 0F,i, 0B,i}: 

tlj{x) = Cn ■ - {n - 

where we have introduced the vectors of the coefhcients 
in the nth unit ceh Cn = {An,Bn,Cn,Dn}- The TM T 
is the matrix describing the change of the amplitudes of 
il) when moving between successive unit cells: 



+ 1 



(4) 



With Q.{x) being the matrix of the four component solu- 
tions, T can be written as [2] 

T = ^{^)-^^{l). (5) 

Since Vt obeys the Dirac equation, it has a constant deter- 
minant, therefore dei[T\ = det[(^(0)"^l^(^)] = 1. Thus, 
the product of eigenvalues of T is one. For a 4 x 4 TM, 
it is useful to specify how the four eigenvalues multiply 
to one. Because det[r^] is constant, we can fix O and I 
in regions without SOIs, where Vt is block diagonal with 
two sub-blocks which are themselves matrices built of two 
component spinors of Ho- As the sub-block matrices are 
also obeying the Dirac equation, their determinants are 
constant as well. Therefore the eigenvalues of T have to 
multiply to one pairwise. Together with the Born-von 
Karman requirement = 1 — the eigenvalues of T 
must be the TVth root of unity — the band condition can 
be written as 



2cos(a>:^) + 2cos(/>:7) = Tr[7l, 



(6) 



where n and n' are the lattice momenta defined by 
jNi with j integer. In order to obtain the r. h. s. of 
this equation we need to know explicitly the matrices of 
spinors outside and inside the SOI regions, I^n and l^so, 
respectively. With these, Eq. (5) can be written as [2] 

r = Fo-^l]N'(0)^N(O^N'(«)^so(a)^^so(0)l^N(0)Fo. 

Here, is a constant 4x4 matrix describing the ampli- 
tudes in the nth unit cell. This matrix does not enter in 
the band conditions but cancels during the calculations. 
In the normal regions the matrix I^n is: 



sgn(£:)e^^^^ — sgn(6:)e ^^^^ 



'l2, (7) 



where (\) = aict8in{ky / kx) in the system without SOI. 
Here, kx^{ky) is the wavevector outside the SOI potential 
in X, (?/)-direction and e the energy of the quasiparticles. 
Finally, I2 is the 2x2 identity matrix in the spin subspace. 
Note, that the system is translational invariant along the 
y direction. In the SOIs regions the matrix l^so reads 



so 



/ie-^^+ 
i(i+ 



V 



ie^^+ 
— i(i+ 
— (i+ 



-le 
—id- 
d- 



-le 
id- 



V. 



'O ■ 



(8) 



In the following, a = ±1 refers to the two spin eigen- 
modes allowed by finite SOIs, and = arctan(/c^/i^cK) 
are their angles of propagation with respect to the x-axis. 
The wavevectors in ^-direction are obtained by eval- 
uating the eigenvalues ofH = Ho + + Hr for uniform 
SOIs and assuming energy conservation and translational 
invariance along the y direction [24]. They are found to 
read 




Ky . 



(9) 
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Further, we hav e in troduced the spin-dependent factors 
da = {s — A)/ ^/K^\ky=o^ and the diagonal matrix Vo = 
diag[Ar+,A/';,A/'_,A/'*] with A/; = e^^-'^ /.y2{l ^ di). 
Calculating the trace exphcitly for ^ = O, the r. h. s. of 
Eq. (6) becomes 

TT[r]=^J'a{K), (10) 

a 

where we have introduced 

J^a{kx) = - (^J-^ + ~^~) sm{kxb) 

-\- 2cos{Kaa) cos{kxb). (11) 

Here, K a and kx are the wavevectors inside and outside 
of the SOI regions, respectively. Note that the func- 
tion does not depend on terms containing the opposite 
spin index. This comes from the fact that, if T were 
unitary, the trace would be the sum of two cosines. If 
= 0, ah matrices building up T are unitary except for 
r^so- In ^so the first two columns correspond to spinors 
associated with eigenvalues a = +1 while the second two 
columns correspond to spinors associated with the eigen- 
value a = —1. Therefore the first two spinors are orthog- 
onal to the second two but only linearly independent to 
each other. Also the second two spinors are only lin- 
early independent. This prevents T from being unitary 
and is the source of the factor in front of the sines in 
Eq. (11). The special form of the four 2x2 subblocks of 
ftso for = O thereby ensures that the sum of the two 
functions J^a{kx) exactly corresponds to the sum of the 
bands associated with the a = mode (plus-band) and 
the a = — 1 mode (minus-band)[25]. Hence, Eq. (6) can 
be written as two independent equations: 

2cos(/^^) =J^^{ka:) ^^2) 
2cos{i<ii) = T-{kx) 

This spin-dependent TM method, which leads to a decou- 
pling of the band equation (6) for the case (/> = O, repre- 
sents the central analytical results of this work. We have 
verified this result by comparing the bands calculated 
from Eq. (12) with those calculated with the method used 
for discussing finite (j) in the last section of this letter [c. f. 
Eq. (18)]. 



IV. COMPARISON WITH THE DISPERSION 
RELATION FOR UNIFORM SOI 

In order to understand what happens for systems with 
periodic SOIs, let us first review what is known about 
graphene with uniform SOIs. There are three interesting 
regions in the A, A-parameter space: (I) A/ A < 1/2, (II) 
A/A > 1/2 and (III) A/A = 1/2. 

An example of the dispersion relations for uniform 
SOIs in these three regions are shown in Fig. 2(a)-(c). 



The energy dispersions associated with a = ^1 {a = —1) 
are depicted as dashed-blue (solid-red) lines. We denote 
the wave vector for the uniform case as k and the energy 
in the uniform case as E and shortly summarize the fea- 
tures for every case. Panel (a) shows that for case (I) 
the spin degeneracy is completely lifted except for /c = O 
where one of the valence and one of the conduction bands 
with opposite spin eigenstates touch. Panel (b) shows 
case (II) in which the intrinsic SOI is dominating over 
the Rashba one. In this case an energy gap opens at 
/c = O where the plus and the minus conduction bands 
are degenerate. In bot h cases the overall behavior shows 
that the energy dispersions acquire a massive nature. In 
panel (c) we show case (III). The energy bands of the 
plus spin state become degenerate at /c = O while for the 
minus spin state there is an energy gap. Importantly, 
for this case, the plus band keeps its linear dispersion, 
contrary to the minus band that still displays massive 
behavior. 



Now let us compare the dispersion relations for uni- 
form SOIs to band structures for modulated SOIs. The 
plus(minus)-bands associated with a = {a = —1) 
are depicted as dashed-blue (solid-red) lines. In the fol- 
lowing we use the Bloch wavevector k instead oi k. In 
region (I), shown in panel (d), there is no gap between 
the valence plus band and the conduction minus band. 
Contrary to the case with uniform SOIs, the two bands 
cross — the valence band is shifted upwards in compari- 
son to the conduction band. Therefore they now intersect 
and there is a degeneracy for two opposite and finite val- 
ues of n. Different is also that gaps — although small 
— open for both bands at the Brillouin zone boundaries 
due to avoided band crossings, as it is usually the case 
when there is a periodic potential. The size of these gaps 
will be estimated later on. In panel (e) we show region 
(II) with modulated SOIs. At = O the degeneracy of 
the two spin polarized conduction bands is lifted because 
of the different mass terms in the uniform SOI case. The 
plus valence band goes through E = ^ and /c = O in the 
uniform case — therefore = O at this point and the 
plus valance band has to go through zero energy at /t: = O 
in the periodic case as well. Again both bands have gaps 
at the Brillouin zone boundary. For region (III), shown 
in panel (f), this is not the case. Due to the linear dis- 
persion of the plus band backscattering is forbidden and 
the band is gapless even at the zone boundaries [3] . This 
can be seen in the band Eq. (12) as well. For A/ A = 1/2 
the prefactor of the sines of the plus band is two and the 
band factorizes to co^^nt) = cos(i^+a + kxh). Also 
is a linear function of e for this case. Concluding this 
comparison we can say that all general features that we 
find in the band structure for periodic SOIs at = O are 
understood from the dispersion relation of uniform SOIs. 
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(a) A = 1 and A = 1/4 (b) A = 1/4 and A = 1/4 (c) A = 1/2 and A = 1/4 




k 



(d) A = 1 and A = 1/4 (e) A = 1/4 and A = 1/4 (f) A = 1/2 and A = 1/4 
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FIG. 2. Panels (a)-(c) show the dispersion relation for uniform SOIs for difFerent values of A and A. Panels (d)-(f) show the 
dispersion relation for periodically modulated SOIs for the same A and A for comparison. The plus-bands (dashed-blue hnes) 
are associated to the a — mode, the minus-bands (sohd-red hnes) to the a — —1 mode. The length of the SOI potential a 
and the length without SOI h are a — h — 7 . 



V. BAND GAPS AT BRILLOUIN ZONE 
BOUNDARY 

As announced before, Eq. (12) enables us to estimate 
the band gaps at the Brillouin zone boundary. Always 
when the expression on the left hand side exceeds mod- 
ulus two there is a gap in the energy spectrum since the 
equation has no real valued solution. In order to calcu- 
late an approximation for the size of the gaps, we follow 
a perturbative approach on the basis of Ref. [26]. First 
we consider A = O and A <C 1. The gap in this case is 
2|[/x=o|7 where Uk=o = — iA/tt is the leading compo- 
nent of the Fourier expansion of the piece-wise constant 
function modulating the intrinsic SOI. For simplicity we 
assume a = b but we obtain a similar solution for a ^ b. 
The band gap for the case without Rashba SOI is thus 
approximately 2A/7r. When A 7^ O we can substitute 
A = A - aA/2 and £ = £ - aA/2 in Eq. (12). Up to first 
order in the expansion parameter bX this becomes 



2 cos{f<ii) =2 cos{Ka(i) cos{eb) 



(13) 



e- A 



sin(i^r+a) ^m{eb) 



aXb 



K a e- A 



e-A K, 



sm{Kaa) cos{be) 0{bX) 



with Ka{s^A) = Ka- The first order term for 6A <C 1 
is rather small and it will never lead to the opening of 
gaps. Actually, when the zeroth order term approaches 
an extremum as a function of e, the first order term ap- 
proaches zero. Thus it is a reasonable approximation to 
neglect the first order term. Therefore — in the pertur- 
bative approximation — the gaps of the bands for A ^ O 



2(A-aA/2) 



-gap 



(14) 



These values coincide with the numerical results. 



VI. DEPENDENCE ON THE 
SPIN-PRECESSION LENGTH 

Finally we want to examine how the band structure 
depends on the periodicity length of the SOIs. As noticed 
in [24] the transmission through a SOI barrier depends on 
the length of the barrier compared to the spin precession 
length and so does the band structure for periodic SOIs. 
First we derive the spin precession length £so, which is 
the length that a carrier needs to cover in order to go back 
to the initial spin state. In the ballistic limit the spin- 
precession length is expressed by -^so = 27r/{K_ — K+). 
A plot of ^so as a function of E is shown in Fig. 3. Note 
that the behavior of ^so is strongly dependent on A but 
the asymptotic value for large E is not. Therefore we 
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FIG. 3. The dependence of ^so on ^ > O (upper panel) and 
E < O (lower panel) for A = 1 and different A is shown. The 
black dashed line refers to the approximation of ^so for large 
E, c.f. Eq. (15). 




FIG. 4. B and stmcture for A = 1 and A = O and 6 = 7 for 
different lengths of the SO potential. Panel (a): a — 2i^l\ — 
i^o and panel (b): a = tt/A = ^^/2. The dashed-blue lines 
correspond to the plus-band and the solid-red lines to the 
minus-band. 



expand in the limit E ^ A, A. This yields a spin 
precession length independent of A and E: 



^SO 



27r 

T 



(15) 



In Fig. 4 the band structure for a = and the band 
stmcture for a = ^go/^ shown in panel (a) and (b) re- 
spectively. In the first case the plus and the minus bands 
lie on top of each other for large energies, where the ap- 
proximation within was calculated is valid. This is 
easily verified by expanding Eq. (12) in the same limit 
up t o first order. One can show that if a = large 
energies only the pre-factor of the sines in Eq. (12) de- 
pends on a. The pre-factor becomes closer and closer to 
one for increasing e. When a = a similar analysis 

shows that the band equation (12) depends on a — in 
this case a is the overall sign — and the two bands are 
shifted with tt in their oscillation. 



VII. FINITE PROPAGATION ANGLE 

After a rather complete discussion of the band struc- 
tures with ^ = O we turn to the case ^ 7^ 0. As explained 



before, for this case Eq. (12) is not valid. Hence, for finite 
angles we evoke Bloch's theorem, which states that 



(16) 



Thus the only change of the spinor after one unit cell is 
the multiplication by a phase factor e^'^^ — therefore 



+ 1 



(17) 



where I4 is the 4x4 identity matrix. Now, equations (4) 
and (17) can be combined in in the following equation 



(T-B)-Cn =0. 



(18) 



which is the band condition for finite 0. Note that in 
absence of SOIs the method summarized by Eq. (18) is 
equivalent to the TM method [2]. But when SOIs are con- 
sidered and a four-dimensional space is needed, it gives all 
conditions on the eigenvalues of T while the TM method 
gives only conditions on the sum of eigenvalues. 

Figure 5 shows the number of solutions of the band 
equation (18) as a function of e and ky for all the allowed 
values of Hi in the first Brillouin zone. As in the previous 
cases, these results can be explained by looking at the 
dispersion relation for uniform SOI. There, a finite value 



of ky — equivalent to a finite value of 



leads to an 



energy gap ai kx = O even if A = 0. Accordingly, in 
panel (a), where the number of allowed bands for the 
modulated case with A = 1 and A = O is shown, a band 
gap growing with ky is displayed. Another feature, which 
is the same for the uniform SOI case and the present 
one, is that if either of A and A are zero, electron-hole 
symmetry is also preserved for finite ky. This can be seen, 
although only in terms of the number of allowed bands, 
in panel (a) and in panel (b), where either of A and A are 
zero and the electron-hole symmetry is preserved, while 
this is not the case in panel (c) and (d), where both, A 
and A, are finite. 



VIII. SUMMARY 

We have studied the band structure of graphene un- 
der a periodically modulated SOI. In this case spin-full 
fermions have to be considered, transforming the orig- 
inal SU (2) structure of the graphene Hamiltonian into 
a SU(2)(g)SU(2) superstructure. For this case the TM 
method — commonly used for calculating band struc- 
tures — is not applicable. Nevertheless, we have pro- 
posed a generalization of the TM method for the spin- 
dependent case at = 0. The new method is possi- 
ble because of the special form of the matrix of spinors 
inside the SO potential. This permits us to calculate 
analytic equations [c.f. Eq. (12)] for the spin-polarized 
bands, which in turn enables us to estimate the gaps due 
to avoided band crossings at the Brillouin zone bound- 
aries [c.f. Eq. (14)] via a perturbative approach. In ad- 
dition, we analyze the dependence of the band structure 
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(a) A = 1 and A = o (b) A = O and A = 1/2 (c) A = 1/4 and A = 1/4 (d) A = 1/2 and A = 1/4 




FIG. 5. The figures show the number of allowed bands as a function of ky and e. Blue corresponds to four allowed bands, 
turquoise to two and white to zero. The SOI potential length is a = 7 and the length without the potential is 6 = 7. 



on the length of the SOI potential compared to the spin- 
precession length [c. f. Fig. 4]. Here we consider the spin 
precession length in the limit of energies much larger than 
the SOI coupling constants, where it is independent of 
the energy [c. f. Eq. (15)]. Furthermore, we investigate 
the band structures for finite <p within an approach based 
on the Bloch's theorem. This leads to a band equation 
[c. f. Eq. (18)], which we solve numerically. We find that 
the modification to the band structure due to finite an- 
gles (/) can be understood by considering the properties 



of the energy dispersion for the case of uniform SOI. 
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